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ABSTRACT 

An  analysis  of  some  ii00,000  station  in^  values  as  reported  in  the  ISC 
bulletin  reveals  substantial  global  variations  in  station  magnitude  bias, 

of  a large 

network  of  stations.  Although  there  are  clear  indications  that  the  biases  are 
functions  of  both  source  region  and  time,  they  appear  to  be  well  correlated 
with  tectonic  structure  and  lateral  variations  in  attenuation  characteristics 
in  the  upper  mantle  under  the  station.  Application  of  these  biases  as  station 
magnitude  corrections  reduces  the  scatter  in  m^  observations  for  a single 
event.  Changes  in  station  distribution  with  time  are  shown  to  introduce 
perceptible  temporal  changes  in  the  shape  of  magnitude-frequency  curve;  these 
can  be  greatly  reduced  by  application  of  the  station  corrections.  These 
corrections,  through  their  effect  on  the  diagram  and  the  magnitude-yield 

relation,  are  applicable  to  the  problems  of  seismic  discrimination. 


defined  as  the  mean  difference  between  station  and  the  average 
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STATION  MAGNITUDE  BIAS  - ITS  DETERMINATION,  CAUSES,  AND  EFFECTS 
I.  INTRODUCTION 

The  existence  of  large  lateral  variations  in  the  attenuation  of  seismic 

waves  has  long  been  recognized.  Evidence  from  the  propagation  of  crustal  body 

12  3 

wave  phases  ’ , the  transmission  of  long  period  teleseismic  P and  S waves  , and 

h 

surface  wave  amplitudes  has  demonstrated  the  existence  of  large  differences  in 
Q in  the  upper  mantle.  These  results  all  indicate  that  attenuation  is  highest 
in  the  regions  of  the  mid-ocean  ridges,  concave  sides  of  island  arcs,  and 
’rift’  structures  such  as  the  western  US,  and  lowest  in  stable  regions  such  as 
shields  and  deep  ocean  basins.  High  attenuation  further  appears  to  be  well 
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correlated  with  high  heat  flow  and  also  certain  negative  velocity  anomalies  ’ . 

In  the  present  work  the  effects  of  these  variations  in  Q upon  body-wave 
magnitude  m^  are  studied.  Substantial  station  biases  have  previously  been 

56  T 

noticed  by  various  authors^’  . Bune  et  al  have  compared  body  wave  magnitudes 
given  by  the  USCGS  (PDE)  and  Russian  sources,  and  found  that  the  PDE  values 
were  substantially  lower.  They  ascribed  this  to  the  regular  contribution  to 
the  PDE  catalog  of  certain  stations  in  the  western  US  which  consistently 
reported  lower  magnitudes  than  most  other  stations. 

The  data  used  here  are  the  station  m^  values  reported  in  the  Bulletin  of 
the  International  Seismological  Centre  (iSC)  for  the  period  196H-73.  Station 
magnitude  biases  are  calculated  for  over  100  stations  and  shown  to  be  well 
correlated  with  the  lateral  variations  in  Q determined  by  previous  authors  by 
other  means.  These  biases  are  then  applied  as  corrections  to  station 
magnitudes  and  this  is  shown  to  achieve  a noticeable  reduction  in  the  scatter  of 
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in^  measurements  for  a single  event.  The  existence  of  such  biases  is  shown  to 
cause  detectable  differences  in  the  shape  of  the  magnitude-frequency  curve, 
temporal  variations  in  the  latter  being  produced  by  changes  in  the  set  of 
stations  reporting  magnitudes.  The  station  corrections  obtained  are  also 
valuable  in  the  context  of  seismic  discrimination  through  their  effect  on  the 
diagram  and  also  in  the  determination  of  magnitude  yield  relationships. 


2 


II.  MAGNITUDE  DATA 


The  ISC  Bulletin  reports  values  of  log  (A/T),  where  A and  T are  the 
amplitude  and  period  of  the  dominant  P-wave  arrival.  Individual  station 
magnitudes  can  then  be  computed  through 

= log  (A/T)  + f(A,h) 

where  f(A,h)  is  a factor  correcting  for  the  source  depth  h and  the  source- 

g 

receiver  separation  A.  The  correction  of  Gutenberg  and  Richter  has  been 
accepted  as  standard  by  most  seismological  organizations,  including  the  NEIS 
(National  Earthquake  Information  Service)  and  ISC,  for  f(A,h).  The  mean  of  all 
reported  determinations  is  then  taken  as  the  event  magnitude.  The  ISC 
bulletin  gives  an  event  m^  only  when  there  are  at  least  3 station  reports: 

the  PDE  does  not  impose  this  restriction. 

9 

Evernden  found  it  necessary  to  adjust  the  distance-depth  correction  of 
Gutenberg  and  Richter  at  distances  of  less  than  20^^:  it  seems  probable  that 

such  adjustments  will  be  highly  region-dependent  and  thus  only  station  reports 
in  the  distance  range  21-100°  have  been  used  here.  With  this  small 
restriction,  the  ISC  bulletins  for  196U-73  contain  UoU,29^  station  log  (A/T) 
reports  for  59^895  events:  of  these  events  Uo,353  had  more  tlian  3 station  log 

(A/T)  reports  and  were  assigned  magnitudes  by  the  ISC.  The  Gutenberg-Richter 
correction  has  been  applied  to  calculate  station  m^  from  the  values  of  log 
(A/T),  and  only  events  satisfying  the  criteria  of  3 or  more  station  reports 
have  been  used  here. 

These  ^40,353  events  had  37^9981  associated  station  magnitude  reports, 
contributed  by  over  500  stations.  Many  of  the  latter  reported  very 
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infrequently  and  are  thus  of  little  use  in  the  bias  calculations.  Only  those 
stations  which  reported  more  than  200  events  in  any  one  year  (a  constraint 
which  requires  only  that  it  report  ^5^  of  seismicity)  were  selected.  There 
were  only  12  such  stations;  Table  I lists  the  number  of  m^  observations  at  each 
during  each  year  of  the  period  196^-73.  The  geographical  location  of  these 
stations  is  shown  in  Figure  1 and  given  in  the  appendix.  The  restriction  of  m^ 
reports  to  these  stations  reduces  the  total  number  of  observations  used  to 
307,^82  for  the  data  period  196^-73.  After  the  station  data  set  has  been 
reduced  to  these  stations  the  event  magnitudes,  defined  as  the  average  of  all 
station  reports,  have  been  recomputed.  The  total  number  of  events  (still 
requiring  3 or  more  observations)  is  reduced  from  UO,353  to  38,3l6.  Figure  2 
illustrates  the  frequency-magnitude  distribution  of  the  events  prior  to  and 
after  the  72-station  restriction:  the  change  in  this  is  clearly  small. 

Surprisingly,  most  events  appear  to  have  been  lost  at  the  higher  magnitude  end 
of  the  distribution;  the  curve  for  the  72-station  magnitudes  is  however  clearly 
smoother  than  that  for  the  original  data.  These  larger  events  which  have 
disappeared  may  be  regarded  as  somewhat  dubious,  since  no  three  of  our  12  best 
stations  have  reported  magnitudes  for  them.  Of  the  U events  of  magnitude  6.7 
on  the  72-station,  3 or  more  station  reports,  event  distribution  only  one  is  an 
earthquake;  the  other  3 are  nuclear  explosions  (Cannikin  and  2 in  Novaya 
Zemlya).  This  may  indicate  an  upper  limit  to  rm  of  ^^^6.7  for  earthquakes. 


Fig.  1.  Location  of  station?  used  in  this  study  (Tables  II  and  IV). 
Also  shown  are  the  boundaries  of  the  seismic  regions  used  in  section 
V,  with  the  total  number  of  events  in  each  reported  by  15  or  more  of 
the  72  stations  of  Table  II. 
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Fig.  2.  Histograms  of  biases  for  BMO  (Oregon)  and  KJN  (Finland)  from  events 
reported  by  10,  15  and  20  stations,  with  the  number  of  events  , mean  p and 
standard  deviation  o’  of  each  distribution. 
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TABLE  I 


Yearly  numbers  of  events  (reported  by  3 or  more  of  72-station  network) 
reported  at  each  of  the  72  stations.  Last  U rows  give  total  of  reports  by 
72  stations,  and  all  stations;  and  events  (requiring  3 or  more  station  m^ 
reports)  for  72-station  and  all-station  networks. 
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TABLE  I (Continued) 


196^4 
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1966 
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3U5 

2U6 
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169 
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99 
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53 
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0 
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0 
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26 
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26J; 

195 
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87 
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35 
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182 
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210 
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5I 
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0 
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0 

0 
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0 

180 

180 

631 
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0 

11 

56 

55 

25 
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265^^ 

3239 

2017 

2127 

2286 

778 
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0 

Total 

reports 


72  sta. 

27019 

35005 

2S811 

29261 

33870 

35268 

323I7 

26832 

26686 

323I7 

All  other 
stations 

2319 

I4I48I 

5112 

8197 

8928 

7377 

7120 

5087 

91 9I 

IO38I 

Events 
(72  st) 

3916 

5051 

3755 

3791 

I023 

I0I5 

3595 

3065 

3232 

381 0 

Events 
(all  St) 

3970 

5132 

3891 

1275 

I503 

I25I 

3696 

3170 

3I66 

3995 

* LAO  (LASA)  was  not  as  bad  as  it  would  appear  to  be  for  1968;  during  this 
year  many  reports  were  assigned  by  the  ISC  to  individual  subarrays 
(e.g.  LFl  etc). 
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III.  BIAS  DETERMINATIONS 


The  bias  b.  , at  the  station  for  the  event  is  calculated  as 
ij 

b. . = m. . - m. 

iJ  ij  J 

where  is  the  station  magnitude  and  m.  the  event  magnitude,  defined  as 


m 


1 r . . 

= — I m ; N = no.  of  stations  reporting. 


J N ij 


Clearly  the  bias  values  will  only  be  significant  if  N is  sufficiently 

large.  A suitable  test  to  find  the  minimum  value  of  N(N  . ) is  to  increase  it 

min 

until  higher  values  cause  no  significant  change  in  the  shape  of  the 

distribution  of  the  biases.  Figure  3 shows  histograms  of  magnitude  biases  for 

various  values  of  N . for  stations  BMO  (Western  U.S.)  and  KJN  (Finland). 

min 

These,  and  similar  diagrams  for  many  other  stations,  show  that  N . = 15  is 

min 

sufficient:  larger  values  reduce  the  size  of  the  data  base  with  neither 

significant  changes  in  mean  value  nor  reductions  in  variance.  The  restriction 
that  N 15  reduces  the  number  of  available  events  to  U668,  with  102,759 
associated  station  m^  reports.  Figure  3 shows  the  distribution  of  these  events 
with  magnitude;  none  are  smaller  than  U.5  and  85^  are  of  m^  ^ 5*0. 

For  each  of  the  72  stations  the  distribution  of  biases  with  respect  to 
event  magnitude  has  been  calculated  for  all  events  reported  (N  ^ 15)  and  the 
histogram  of  the  biases  plotted.  Histograms  for  9 stations  are  shown  in  Figure 
U.  For  all  the  stations,  the  normal  distribution  is  a remarkably  good 
approximation  to  that  observed,  and  thus  only  the  mean  and  its  associated 
standard  deviation  are  required  to  characterize  the  nature  of  the  bias 
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10000 


Fig.  3.  Incremental  magnitude-frequency  curves  (each  point  at  a given 
magnitude  is  number  of  events  of  magnitude  m^ 

(a)  all  events  with  m^  reported  by  any  station 

(b)  all  events  reported  by  3 or  more  stations 

(c)  all  events  reported  by  3 or  more  of  the  72 
stations  chosen 

(d)  all  events  reported  by  15  or  more  of  the  72  stations  chosen- 
Time  period  is  196^-73- 
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NUMBER  EVENTS 


Fig.  4.  Histogram  of  biases  for  9 stations,  data  base  196U-73- 
Also  given  are  the  number  of  events  in  the  sample  and  the  mean 
and  standard  deviation  of  the  distribution. 
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distribution.  Table  II  gives  the  number  of  events  reported  (requiring  N i 15) 
at  each  station  and  corresponding  values  of  mean  bias  b,  standard  deviation 
and  error.  The  error  as  given  here  is  the  standard  deviation  of  the  mean, 
defined  as 


1 


b 


15 

I (b.-b)^ 
J=1 


being  the  number  of  events  with  more  than  15  associated  station  ni^ 
values  reported  at  the  station. 

The  mean  biases  calculated  can  be  seen  to  range  from  -0.32  to  +0.37  m^ 
units  and  the  standard  deviations  between  0.2^  and  0.50;  in  only  8 cases  (BNG, 
CAN,  EUR,  HFS,  LAO,  TSK,  WIN)  does  it  exceed  O.^i.  The  mean  itself  is  well 
estimated  because  of  the  large  sample  size:  the  error,  a weighting  of  the 

standard  deviation  by  the  sample  size,  is  a measure  of  the  accuracy  of  the 
mean . 

Possible  causes  of  these  large  standard  deviations  include  temporal 
changes  in  station  bias  characteristics,  severe  dependence  upon  source  region 
due  to  extreme  lateral  variations  in  either  structure  or  attenuation  beneath 
the  station  or  source,  and  a possible  dependence  of  bias  upon  event  magnitude 
through  either  detection  characteristics  or  changes  in  reported  dominant 
period,  and  thus  attenuation,  through  source  properties.  Each  of  these  factors 
is  discussed  below. 
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TABLE  II 


Number  of  events  reported  at  station  of  T2-station  network,  requiring  3 or 
more  {N3)  and  15  or  more  (N^^^)  reports  per  event.  Columns  3,  and  5 give 
mean  (b),  standard  deviation  (o^)  of  distribution  of  biases  for  events 

and  error  (defined  in  text  as  in  estimation  of  mean  biases.  Last 

column  gives  structural  type  assigned  to  each  station  site  (S-Shield,  P- 
Platform,  R-Rift,  O-Oceanic,  and  F-Foldbelt  and  Seismic). 


"15 

b 

% 

error 

X 10^ 

Structure 

ALE 

1390 

1160 

-o.oi* 

0.29 

0.85 

S 

ALQ 

5205 

1917 

-0.20 

0.33 

0.77 

R 

ASP 

1813 

1*66 

-0.05 

0.35 

1.59 

S 

BHA 

2829 

980 

-0.28 

0.32 

1.03 

R 

BKS 

i1)T6 

877 

+0.18 

0.38 

1.31 

P 

BMO 

26192 

3777 

-0.29 

0.35 

0.  56 

R 

BNG 

39l<U 

1029 

-0.07 

0.50 

1.56 

P 

BNS 

1975 

1695 

+0.20 

0.29 

0.71 

P 

BOZ 

27l<8 

978 

-0.06 

0.31 

1.00 

R 

BUL 

*^793 

1368 

-0.07 

0.29 

0.78 

R 

CAN 

21j62 

507 

-0.02 

0.1*0 

I.7I* 

F 

CAR 

1289 

67 1* 

+0.13 

0.38 

1.1*6 

F 

CIR 

2623 

938 

-0.27 

0.30 

0.97 

R 

CLL 

2726 

1715 

+0.20 

0.32 

0.78 

P 

CLK 

2lt51 

902 

-0.27 

0.28 

0.93 

R 

COL 

IIU18 

2579 

+0.01 

0.33 

0.65 

F 

COP 

H*57 

1290 

+0.36 

0.26 

0.72 

P 

CPO 

11569 

2510 

-0.07 

0.35 

0.70 

P 

DUG 

6117 

1753 

-0.15 

0.35 

0.83 

R 

EDM 

1709 

1181 

+0.37 

0.28 

0.82 

P 

FKA 

309^ 

1395 

+0.  00 

0.33 

0.89 

P 

EUR 

1U997 

2913 

-O.2I* 

0. 1*0 

0.75 

R 

FUR 

2I98 

1816 

+0.10 

0.38 

0.90 

F 

GDH 

11*39 

1257 

+0.00 

0.31* 

0.97 

S 

GIL 

6138 

181*2 

-o.oi* 

0.35 

0.81 

F 

GOL 

1*058 

1770 

-0.28 

0.39 

0.93 

R 

GRF 

1U5O 

liol* 

+O.2I* 

0.28 

0.85 

P 

HFS 

2917 

895 

+0.05 

0.1*5 

1.50 

S 

HYB 

2106 

1166 

+0.19 

0.37 

1.09 

S 

KEY 

3678 

2218 

+0.02 

0.27 

0.57 

S 

KHC 

1*092 

2720 

+0.10 

0.26 

0.  50 

F 

KJF 

2137 

1036 

+0.09 

0.28 

0.88 

S 

KJN 

691*0 

2l*00 

+0.1I* 

0.30 

0.61 

S 

KOD 

2291* 

1210 

+0.06 

0.31 

0.88 

S 

KON 

Il4l42 

1093 

+0.07 

0.30 

0.91 

S 

KRA 

171*3 

1251 

+0.22 

0.29 

0.83 

P 

KRR 

3159 

IOU7 

-O.2I* 

0.30 

0.9I* 

R 

KTG 

2221 

1798 

+0.02 

0.30 

0.71 

S 

LAO 

10202 

1705 

-0.10 

0.1*7 

1.15 

P 
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TABLE  II  (Continued) 


“3 

^15 

L 

"b 

error 

X 10^ 

LJU 

2501 

1739 

+0.29 

0. 30 

0.73 

LON 

2093 

1333 

-0. 30 

0.37 

1.03 

LOR 

1879 

11)47 

+ 0.06 

0. 12 

I.2I 

LPS 

679 

+ 0.  0^4 

0. 31 

1.31 

MBC 

2327 

1256 

+ 0.1^4 

0.31 

0.97 

MOX 

I46U5 

2762 

+0.02 

0.27 

0. 52 

MUN 

1738 

397 

+0.15 

0. 37 

1.85 

NAO 

3399 

755 

-0.09 

0.29 

1. 07 

NDI 

1653 

906 

+0.33 

0.37 

1.23 

NEW 

2116 

1221 

+0.05 

0.30 

0.86 

NIE 

1566 

1096 

-0.02 

0.  33 

1.00 

NOR 

3693 

1986 

-O.lU 

0.33 

0.75 

NP- 

393)4 

1035 

-0.00 

0.  38 

1.19 

NUR 

7256 

31)49 

+0.19 

0.  30 

0.  5)4 

PMG 

Ul(35 

1116 

0 

1 — 1 

d 

+ 

0.38 

1.12 

PMR 

6016 

2075 

-0.08 

0.37 

0.82 

PNT 

1178 

663 

+0.13 

0.  30 

1.15 

POO 

28I1 

1288 

+0.IT 

0.  36 

1.00 

PRE 

19)47 

802 

-O.OT 

0.  39 

1.39 

PRU 

2779 

2210 

+0.0U 

0.2)4 

0. 51 

RES 

l)t35 

106)4 

+0.13 

0.37 

1.16 

RIV 

1730 

507 

+0.  31 

0.33 

1.50 

SHL 

2753 

769 

+0.11 

0.33 

1.22 

SJG 

15)47 

718 

+0.2U 

0.  38 

1.  )40 

STU 

1752 

1)43)4 

+0.29 

0.31 

0.81 

TFO 

185)46 

2395 

-0.32 

0.35 

0.71 

TRN 

1337 

70I 

+0.07 

0. 35 

1.35 

TSK 

2762 

10I7 

-0.07 

0.I5 

i.lo 

TUC 

3033 

1263 

-0.II4 

0.25 

0.71 

TIJL 

2820 

1152 

+0.21 

0.32 

O.9I 

UBO 

23157 

2828 

( — 1 
1 — 1 

d 

1 

0.  38 

0.72 

WIN 

1259 

6)43 

-0.09 

O.I43 

1.72 

WMO 

13101 

1658 

-0.17 

0.  31 

0.76 

Structure 


F 

R 


P 

S 

s 

s 

R 

F 

S 


F 

F 

F 

S 

S 

F 

S 


P 

F 

F 

P 

R 

F 

F 

R 

P 

R 

S 

P 


ih 
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IV.  TEMPORAL  VARIATION  IN  STATION  BIAS 


Bias  distributions  have  been  calculated  for  each  individual  year  of  data 
as  well  as  for  the  entire  time  period  196^1-73 • In  general  the  mean  station 
biases  show  little  variation  from  year  to  year,  but  for  some  stations,  and 
particularly  for  some  arrays,  dramatic  changes  in  the  bias  distribution  took 
place  with  time. 

Biases  for  any  individual  year  are  calculated  only  when  there  are  more 
than  100  observations  of  events  with  15  or  more  station  m^  reports.  For  12 
stations  the  largest  difference  between  the  mean  bias  for  any  one  particular 
year  and  that  for  the  entire  time  period  196^1-73  exceeded  0.2  m^  units.  These 
stations  are  in  many  cases  the  same  as  those  with  large  standard  deviations  of 
bias  distribution  as  given  in  Table  II.  In  Figure  5 the  change  in  bias 
distribution  is  shown  for  these  12  stations  as  a function  of  date.  The  dots 
indicate  the  mean  bias  and  the  bars  ± one  standard  deviation  of  the 
distribution.  Note  that  the  latter  do  not  indicate  errors  in  the  mean;  the 
mean  is  extremely  well  estimated  because  of  the  large  sample  sizes  (at  least 
100). 

Particularly  alarming  cases  are  those  of  HFS,  LAO  and  EUR.  For  the  first 
two  (arrays  in  Sweden  and  Montana)  there  appears  to  have  been  a severe 
degradation  in  performance,  characterised  by  large  standard  deviations,  in 
certain  years,  and  at  EUR  (Nevada)  there  is  a remarkable  decrease  in  bias  over 
the  ten-year  time  interval.  It  is  difficult  to  conceive  of  any  rational 
explanation  for  the  latter:  nearby  stations  DUG  and  UBO  do  not  exhibit  any 
trend. 
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Fig.  5*  Annual  variations  in  mean  “bias  for  the  12  stations  for 
which  the  mean  in  any  one  year  differs  from  that  for  the  entire 
time  period  (196^-73)  by  more  than  0.2  m-^  units.  The  bar  de- 
notes t one  standard  deviation  cr-j^  of  the  bias  distribution, 
not  the  mean. 
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For  stations  COL  and  LON  the  largest  differences  from  the  average  hias 
occur  in  years  196^^  and  65  - during  these  seismicity  was  greatly  dominated  hy 
events  in  Alaska  (6^4)  and  the  Rat  Islands  (65);  some  regional  source  effects 
may  be  operating  here.  PRE  and  WIN,  both  in  southern  Africa,  show  a sudden 
decrease  in  bias  in  1972  - this  is  not  noticeable  on  similar  plots  for  stations 
in  East  Africa. 

Most  of  these  temporal  changes  defy  any  rational  explanation.  In  general 
the  years  with  largest  deviations  of  mean  bias  are  also  those  with  the  largest 
variance:  thus  the  statistical  significance  of  these  trends  may  be  small.  The 

vast  majority  of  stations  used  here  show  no  such  trends,  including  those  which 
have  reported  most  events  over  the  entire  time  interval  (e.g.,  UBO,  KJN, 

TFO) . 
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V.  SOURCE  REGION  VARIATIONS  IN  BIAS 


We  may  consider  the  mean  station  bias  b as  defined  in  section  III  to 
consist  of  the  following  factors 


b = b" 


cs  us 
b + b 


+ b 


M 


c s c r 

where  b , b are  introduced  by  crustal  structure  at  the  source  and  receiver 
respectively,  b^^  and  b^^  by  upper  mantle  structure  in  the  same  regions,  b^  in 

g 

the  lower  mantle  part  of  the  ray  path,  and  b by  the  source  radiation  pattern. 
All  of  these  factors  clearly  may  depend  upon  the  source-receiver  configuration 
through  source  take-off  angle  and  receiver  arrival  angle.  The  Gutenberg- 

Richter  distance-depth  correction  may  be  considered  to  approximate  the  effects 

cs  , us  , ur  .CTx  ^ ^ ^ ^ 

lb  +b  +b+b  +bj  globally;  the  biases  measured  then  really  measure 

s 

deviations  from  this  average  behavior.  The  term  b accounts  for  deviations 

from  the  average  amplitude  in  a small  region  surrounding  the  source  due  to 

radiation  pattern  effects.  Station  site  effects,  e.g.,  seismometer-ground 

coupling  variations  due  to  the  medium  (hard  rock,  alluvium)  upon  which  the 

cr 

station  is  situated,  are  included  in  b 

Seismic  sources,  by  their  very  nature,  have  a tendency  to  be  located  in 
regions  of  high  lateral  inhomogeneity,  and  thus  it  may  be  expected  that 
features  such  as  the  anomalously  high  attenuation  on  the  concave  side  of  island 
arcs  can  seriously  effect  determinations.  Studies  leading  to  the 
development  of  plate  tectonic  theory  have  indicated  consistency  of  fault  plane 
solutions  and  thus  presumably  radiation  patterns  over  large  source  regions. 
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TABLE  III 

Deviation  of  mean  biases  for  each  source  region  from  mean  bias 
for  all  regions  (column  3,  Table  II),  for  sets  of  contiguous 
stations  in  Germany,  East  Africa,  and  the  western  US. 
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The  events  used  here  have  been  separated  into  11  major  source  regions  and 
biases  computed  in  the  same  manner  as  before  for  events  from  different  source 
areas  recorded  at  a given  station.  As  in  the  pi'evious  section,  100  or  more 
events  are  required  for  determination  of  mean  bias.  Variations  of  bias  with 
source  region  are  generally  somewhat  higher  than  those  with  time  given  in  the 
previous  section.  For  l6  stations  the  mean  bias  determined  from  events  in  a 
given  source  region  differed  by  more  than  0.2  m^  units  from  that  for  all 
sources  given  in  Table  II.  For  many  stations  there  are  an  insufficient  number 
of  events  for  individual  regional  variations  in  biases  to  be  calculated,  and 
thus  comparison  of  biases  from  a number  of  regions  is  severely  limited  to  the 
better  stations.  There  was  no  indication  that  any  one  particular  source  region 
was  more  anomalous  in  terms  of  bias  characteristics  than  any  other. 

Many  of  the  stations  considered  here  happen  to  be  concentrated  in  regions 
which  are  small  in  extent  compared  to  source-receiver  distances  and  we  may 
therefore  expect  to  see  consistencies  in  the  variations  of  mean  bias  with 
source  region.  Three  particularly  small  receiver  regions  are  those  containing 
stations  in  Germany,  East  Africa,  and  the  Western  U.S.  Table  III  gives  the 
variations  of  bias  with  source  region  for  each  of  these  area. 

It  can  be  seen  that,  even  over  these  small  receiver  regions,  there  is 
remarkably  little  consistency  in  variations  of  mean  bias  with  source  region, 
the  only  possible  exception  being  that  for  events  in  region  U (Japan)  to  the 
western  U.S.  Similar  tables  for  other  receiver  areas  such  as  Scandinavia  and 
India  also  fail  to  reveal  any  correlation  of  bias  with  source  region.  Had  we 
selected  smaller  source  regions  it  is  likely  both  that  regional  variations  in 
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bias  and  their  consistency  across  receiver  areas  would  have  been  more 

pronounced,  particularly  in  the  case  of  ray  paths  travelling  down  descending 

lithospheric  slabs:  unfortunately  the  data  base  is  unsufficient  to  test  this. 

This  lack  of  correlation  may  be  taken  to  indicate  that  the  tei'ms 
c s u.  s M 

(b  *’  + b + b ),  however  large  they  may  be  for  individual  ray  paths,  tend  to 
average  out  in  such  a manner  that  they  cannot  be  resolved  in  the  present  study. 
It  also  indicates  that  the  Gutenberg-Richter  distance-depth  correction  is  not 

grossly  in  error.  In  particular,  much  of  the  observed  station  bias  may  be  due 

ur  . cr 

to  the  term  b , and  its  variation  with  source  region  to  b . This  does  not  of 

c s us  c r ur 

course  imply  that  b and  b are  not  as  large  as  b and  b : in  fact  the 

grevater  lateral  heterogeneity  in  source  regions  probably  means  that  they  will 

be  larger;  but  they  cannot  be  resolved  by  the  present  means. 
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VI.  VARIATION  OF  BIAS  WITH  MAGNITUDE 


There  exist  two  plausible  reasons  for  expecting  a variation  of  station 
bias  with  magnitude.  As  mentioned  in  section  II,  these  are  the  variation  of 
the  frequency  content  of  seismically  radiated  waves  with  source  size,  and  the 
station  detection  characteristics. 

All  present  theories  of  the  seismic  source  incorporate,  in  various 
manners,  an  increasing  proportion  of  energy  at  longer  periods  as  the  source 
size  increases.  This  therefore  implies  that  not  only  the  amplitude,  but  also 
the  dominant  period,  of  the  initial  P-wave  arrival  from  which  is  measured, 
increases  with  the  size  of  the  seismic  source.  If  the  attenuation  of  seismic 
waves  was  laterally  homogeneous,  then  its  dependence  upon  frequency  could  not 
be  detectable  by  a study  of  m^  biases,  since  the  average  itl^  with  respect  to 
which  we  measure  biases  would  depend  upon  frequency  in  the  same  manner.  Since 
attenuation  is  clearly  laterally  dependent,  as  shown  by  the  biases  in  Table  II 
and  the  other  studies  mentioned  in  the  introduction,  we  may  expect  a variation 
of  bias  with  magnitude  through  the  dependence  of  attenuation  as  a function  of 
period  for  average  Q in  a given  region.  Unfortunately,  we  cannot  directly 
measure  the  variation  of  bias  with  period  since  the  ISC  Bulletins  generally 
give  log  (a/T)  and  not  A and  T individually.  We  may,  however,  be  able  to 
detect  some  variation  of  bias  with  event  magnitude. 

Ringdal^^,  Christoffersson  et  al.^^  and  others  have  considered  the 
effects  of  station  detection  thresholds  upon  magnitude  determinations.  Their 
models  incorporate,  in  a complicated  manner,  the  effects  of  both  station  biases 
and  detection  thresholds  upon  the  relationship  between  station  itl  and  ^true* 
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The  joint  estimation  of  the  relevant  parameters  for  each  station  is  an 


extremely  involved  procedure.  Figure  3 shows  the  magnitude-frequency 
distribution  of  the  events  used  in  the  bias  calculations.  Ringdal^^  has 
estimated  the  detection  capability  of  the  Norwegian  seismic  array  (NORSAR)  for 
events  in  the  Japan-Kuriles-Kamchatka  region  and  finds  a 90^  detection 
capability  at  m^  = ^.27-  NORSAR  is  one  of  the  better  stations  studied  here, 
and  it  seems  certain  that  90^  detection  capabilities  for  many  of  the  stations 
will  be  higher.  Figure  3 shows  that  just  over  85/?  of  the  events  used  in  the 
bias  calculations  have  magnitudes  > 5«0*,  this  may  still  be  below  the  90% 
detection  capability  of  some  stations,  particularly  since  biases  are  calculated 
for  all  stations  which  report  values  for  200  or  more  events  annually.  Many 
of  the  stations  in  Table  I have  reported  no  more  than  10%  of  the  total  number 
of  events  in  Figure  3.  The  percentage  of  total  events  reported  is  however 
clearly  a function  of  the  geographical  distribution  of  seismicity  with  respect 
to  the  station:  for  example  stations  in  the  Western  U.S.  can  observe  most  of 

the  circumpacif ic  seismicity  - those  in  Europe  and  Africa  are  not  so  fortunate 
in  this  respect. 

In  order  to  ascertain  whether  these  effects  are  serious,  the  data  base 
used  has  been  separated  into  3 magnitude  classes:  ^.5  ^ m^  < 5-0, 

5-0  < < 5-5 

station  for  events  in  each  magnitude  class,  requiring  once  again  100  or  more 
measurements  to  determine  a mean  bias.  The  variation  of  mean  bias  with 
magnitude  is  small:  in  only  13  cases  does  it  exceed  0.1  iri^  units  and  even 

these  are  of  dubious  statistical  significance.  It  may  reasonably  be  assumed 


and  5 • 5 < 


< D.O.  Biases  have  heen  calculated  for  each 
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that  one  measure  of  detection  capability  is  the  number  of  events  reported  by  a 
particular  station,  but  there  appears  to  be  no  significant  correlation  between 
this  and  the  variation  of  mean  bias  with  magnitude. 

The  restriction  of  bias  calculation  to  events  with  15  or  more  station 
reports  requires  that  the  sources  used  are  of  fairly  large  size  and  at  this 
level  we  may  be  above  the  detection  thresholds  of  most  stations.  The  effects 
of  such  thresholds  is  probably  serious  at  smaller  magnitudes,  for  which  the 
biases  obtained  may  not  be  valid:  the  only  measure  of  this  is  whether 

application  of  the  biases  at  these  levels  reduces  the  scatter  in  magnitude 
observations  and  thus  the  variance  in  average  magnitude.  This  will  be  exfimined 
in  a later  section. 
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VII.  REGIONAL  VARIATIONS  IN  BIAS 


In  the  previous  three  sections  the  dependence  of  station  magnitude  bias 
upon  time,  source  region,  and  magnitude  has  been  examined.  In  some  cases  the 
variation  of  mean  bias  with  time  has  been  shovn  to  be  associated  with  large 
variances  of  the  distribution,  and  the  effects  of  magnitude  variations  appear 
to  be  small.  It  has  been  tentatively  concluded  in  section  V that  the  mean 
biases  obtained  do  represent  the  effects  of  attenuation  in  the  region  near  the 
receiver.  In  this  section  we  examine  the  correlation  between  these  biases  and 
other  indications  of  lateral  variations  in  attenuation. 

Figures  6 through  12  indicate  the  variations  in  mean  bias  (hereafter 

referred  to  as  station  bias)  across  North  America,  Europe,  Africa,  Australia, 

and  India.  In  the  continental  U.S.  (Figure  T)  the  large  differences  in 

attenuation  between  the  western  and  eastern  U.S. , previously  noted  by  Romney 
1 3 ^ 

et  al.  , Solomon  and  Toksoz  and  Solomon  , are  clearly  apparent.  Attenuation 
is  obviously  higher  in  the  western  U.S.  than  in  the  older  stable  regions  of  the 
east . 

Unfortunately  our  original  data  base  of  72  stations  contains  only  3 in 
the  eastern  U.S. A.,  and  to  improve  this  situation  biases  have  been  calculated 
for  30  further  stations  including  some  in  this  area.  These  stations  do  not 
satisfy  the  criteria  (viz  at  least  200  reports/year)  of  section  II;  however 
each  reported  at  least  500  observations  over  the  time  interval  196^-73. 

Because  of  the  paucity  of  data,  we  are  unable  to  determine  any  temporal  or 
source  region  variations  in  these  biases,  and  they  must  necessarily  be 
considered  less  reliable.  Table  IV  lists  these  stations,  the  number  of 
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Fig.  6.  Mean  biases  for  stations  in  the  continental  USA. 
Values  in  parentheses  are  from  Table  IV;  all  others  are 
from  Table  II. 
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Fig.  7.  Mean  biases  for  stations  in  Canada,  Alaska,  and  Greenland, 
in  parentheses  are  from  Table  IV;  all  others  are  from  Table  II. 


Values 
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Fig.  8.  Mean  biases  for  stations  in  Europe.  Values  in  paren- 
theses are  from  Table  IV;  all  others  are  from  Table  II. 
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Fig.  9-  Mean  biases  for  stations  in  East 
parentheses  are  from  Table  IV;  all  others 


Africa.  Values  in 
are  from  Table  II. 
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Fig.  10.  Mean  biases 
theses  are  from  Table 


for  stations  in  India.  Values  in  paren- 
IV;  all  others  are  from  Table  II. 
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Fig.  11.  Mean  biases  for  stations  in  Australia.  Values  in 
parentheses  are  from  Table  IV;  all  others  are  from  Table  II. 
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Fig.  12.  Mean  biases  for  stations  in  French  Polynesia.  Large 
island  at  lower  left  is  Tahiti.  Values  in  parentheses  are  from 
Table  IV;  all  others  are  from  Table  II. 
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TABLE  IV 


Number  of 

events  1 

[assigned  m^^ 

values  by  15  or 

more  stations) 

reported 

by  30 

supplementary  stations. 

Next  columns  give  mean  and  standard 

deviation 

of  bias 

distribution 

i for  these  events,  error,  and 

structural 

type 

(for 

definitions  of  last 

two  see  caption 

to  Table  II ) . 

N(15) 

b 

a. 

error 

Structure 

b 

X 10^ 

ABU 

702 

+0.29 

0.53 

2.0I* 

F 

ATL 

213 

+0.35 

0.51 

3.1*0 

P 

BRG 

7U9 

-0.11 

0.26 

0.96 

F 

BUH 

klh 

-o.oU 

0.  38 

1.90 

F 

DAG 

39^ 

-0.02 

0.27 

1.35 

F 

ESK 

650 

+0.19 

0.37 

1.1*2 

P 

FFC 

568 

+0.08 

0.29 

1.21 

S 

FLO 

k95 

+0.26 

0.  31 

1.1*1 

P 

GBA 

350 

+ O.OI4 

0.37 

1.95 

S 

ILG 

31 

+O.OI4 

O.U2 

2.1*7 

S 

INK 

90U 

+0.15 

0.29 

0.97 

S 

KBL 

578 

+0.09 

0.31 

1.29 

F 

KLG 

639 

-0.02 

0.1*5 

1.80 

S 

LHN 

hj6 

+0.13 

0.36 

1.61* 

S 

LPB 

+0.07 

0.31 

1.1*8 

F 

MAW 

190 

+0.11 

0.35 

2.50 

S 

MTD 

h32 

-0.22 

0.  30 

1.1*3 

R 

OIS 

392 

-0.08 

0.1*1 

2.13 

F 

PAE 

2U8 

-0.03 

0.29 

1.81 

0 

PMO 

277 

+0.08 

0.35 

2.06 

0 

PNS 

370 

-0.08 

0.1*5 

2.37 

F 

PPN 

233 

-0.19 

0.  32 

2.13 

0 

PPT 

237 

+0.06 

0.29 

1.93 

0 

RAB 

766 

+0.12 

0.1*3 

I.5I* 

F 

RCD 

26k 

+O.2U 

0.31 

I.9I* 

P 

ROL 

392 

+0.19 

0.31 

1.55 

P 

RUV 

227 

+0.09 

0.37 

2.1*6 

0 

TPT 

299 

+0.06 

0.  3I* 

2.00 

0 

TVO 

235 

+0.08 

0.28 

1.87 

0 

VAH 

283 

+0.02 

0. 32 

1.88 

0 
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observations  of  events  with  at  least  15  associated  reports  given  by  each, 
and  the  corresponding  mean,  standard  deviation  and  error.  Biases  for  these 
stations  are  included  in  Figures  7 through  12;  they  are  given  in  parentheses. 

Inclusion  of  stations  ATL,  FLO,  ROD  and  ROL  further  clarifies  the 
differences  in  attenuation  in  the  U.S.  We  have  only  one  station  in  California 
(BKS);  its  positive  station  bias  does  agree,  however,  with  a decrease  in 
attenuation  near  the  west  coast  noted  by  other  authors.  Differences  in  station 
bias  of  up  to  0.6  ra^  units  are  apparent  in  the  U.S. 

An  examination  of  Figures  7 through  12  confirms  that  station  biases  are 

highest,  and  thus  attenuation  lowest,  in  shield  regions  such  as  Canada,  India, 

Scandinavia  and  Australia.  The  only  region  where  biases  as  low  as  those  in  the 

western  U.S.  are  observed  is  East  Africa  (Figure  9):  the  effect  of  the  East 

12  13 

African  rift  valley  is  apparent  and  surface  wave  dispersion  studies  ’ have 
shown  clear  similarities  in  velocity  structure  between  these  two  regions. 

We  have  8 stations  in  an  oceanic  region:  these  are  all  located  in  the 

small  area  in  French  Polynesia  shown  in  Figure  12.  Seven  of  these  stations 
have  mean  biases  in  the  range  -.03  to  +.09  m^  units.  Other  small  regions  shown 
in  Figures  6 through  11  also  show  a consistency  in  mean  bias  (Rhodesia,  Figure 
9;  Northern  Germany  and  Finland,  both  Figure  8). 

Each  station  in  Tables  II  and  IV  has  been  assigned  to  one  of  5 tectonic 
structures:  shield,  aseismic  platform,  rift,  oceanic,  and  foldbelt  (including 

present  seismic  regions).  These  are  denoted  by  S,  P,  R,  0,  and  F respectively 
in  Tables  II  and  IV.  Figure  13  shows  histograms  of  mean  station  bias  for  each 
region.  The  distinction  between  rift  structure  and  all  the  other  types  is 
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Fig.  13. 
classes : 
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Mean  biases  separated  into  sub-station  structure 
structure  types  from  Tables  II  and  IV. 
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clear;  that  between  shield,  platform,  foldbelt  and  oceanic  is  not  so  obvious. 

The  sharp  peak  in  the  oceanic  biases  is  certainly  not  representative;  all  these 

biases  are  from  a very  small  region  in  the  South  Pacific  (see  P'igure  12). 

Biases  in  platform  regions  are  slightly  higher  than  in  shield:  tliis  is  not  too 

consistent  with  other  measurements  of  attenuation  for  these  two  structural 

types.  In  seismic  regions,  particularly  on  the  concave  side  of  island  arcs, 

2 

Molnar  and  Oliver  have  postulated  high  attenuation  on  the  basis  of 
propagation  characteristics;  the  histogram  of  mean  biases  for  stations  in 
seismic  regions  (mainly  from  Japan,  South  America,  Alaska  and  the  Caribbean)  is 
not  consistent  with  this.  A comparison  of  the  geographical  distribution  of 
seismic  activity  with  that  of  these  stations  reveals  that  the  region  directly 
behind  and  above  the  descending  lithosphere  in  their  vicinity  is  poorly  sampled 
by  teleseismic  ray  paths:  this  may  be  a contributing  factor  to  these 

unexpectedly  high  (or,  more  precisely,  non-low)  values  of  station  bias. 

Despite  all  these  reservations,  the  correlation  between  tectonic  type  and 
bias  shown  here  is  sufficiently  good  that  our  earlier  conclusion  that  the 
biases  as  measured  reflect  upper  mantle  conditions  near  the  receiver  would 
appear  to  be  justified.  The  agreement  shown  is  not  unexpected  in  view  of 
previous  studies,  as  mentioned  in  the  introduction;  it  is  however  gratifying 
that  such  a poor  measure  of  amplitude  as  can  reveal  some  of  these 
differences  in  attenuation. 
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VIII.  CORRELATION  OF  BIASES  WITH  VELOCITY  ANOMALIES 


An  association  of  regions  of  high  attenuation  with  negative  velocity 

1 5 it 

anomalies  has  been  remarked  by  many  previous  authors  ’ ’ . Surface  wave  phase 

velocities,  as  summarised  by  Knopoff^^  are  highest  in  shield  and  aseismic 

platform  regions  and  lowest  in  rift  and  foldbelt  areas,  these  being  correlated 

with  high  and  lower  upper  mantle  velocities  respectively.  Marshall^^  has 

measured  biases  for  many  of  the  stations  used  here:  these  agree  substantially 

with  those  given  in  Table  II.  He  has  demonstrated  a relation  between  station 

bias  (and  thus  attenuation)  and  P velocities  beneath  the  station  for  the 

n 

continental  U.S. , and  assumed  this  relation  to  hold  elsewhere  in  the  world. 

Figure  lUa  shows  magnitude  biases  measured  for  stations  within  the  U.S.  versus 

P velocities  beneath  the  stations.  The  P values  used  here  are  from  the  map 
n n 

17 

of  Herrin  . There  is  no  doubt  that  bias  increases  with  increasing  P^, 

though  there  are  some  anomalies  (particularly  V/MO,  about  which  more  is  said  in 

1 6 

section  X).  Marshall  has  applied  the  P^-bias  relation  he  derives  to  improve 

the  magnitude-yield  curve  for  explosions  and  demonstrated  that  the  use  of  both 

receiver  and  source  biases  can  dramatically  improve  the  linearity  of  this 

curve.  Figure  lUb  shows  biases  versus  P travel-time  station  anomalies  for  the 

1 8 

U.S. , from  Figure  20a  of  Hales  and  Herrin  . Here  again  the  correlation  is 
almost  convincing,  if  we  ignore  CPO  and  WMO  (for  possible  justification  see 
section  X). 
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Fig.  ih , Mean  biases  for  continental  US  stations  versus 

(a)  sub-station  velocity  (from  Reference  IT) 

(b)  station  P travel-time  anomalies  (from 
Ref erence  l8 ) . 


IX.  APPLICATION  OF  BIASES  AS  STATION  CORRECTIONS 

As  veil  as  their  intrinsic  geophysical  interest,  a clear  application  of 
the  biases  ve  have  obtained  is  their  use  as  station  corrections  to  reduce 
scatter  in  the  determination  of  m^  for  a particular  event.  We  may  correct  each 
individual  station  ni^  by 

- bias 

This  method  is  analogous  to  the  joint  epicentre  determination  technique  of 
Douglas^^  vhich  uses  travel  time  residuals  from  a large,  veil-recorded  (master) 
event  as  station  travel-time  corrections  for  smaller  events  nearby.  The 
epicentres  are  then  relocated  relative  to  the  master  event.  The  average 
(event)  iri^  values  ve  obtain  vill  also  be  corrected  so  that  they  are  more 
relatively  accurate;  their  absolute  values  are  undetermined. 

t 

Ideally  ve  should  also  apply  a bias  correction  for  attenuation  in  the 

vicinity  of  the  source,  since,  as  shovn  in  sections  V and  VII,  the  station 

biases  as  measured  here  reflect  mainly  attenuation  in  the  vicinity  of  the 

station.  Most  earthquakes  occur  in  subduction  zones  and  here,  as  the  evidence 

2 

of,  amongst  others,  Molnar  and  Oliver  has  shovn,  attenuation  is  high  on  the 
concave  side  of  the  Pacific  island  arcs.  Reciprocity  implies  that  station  bias 
in  seismic  regions  can  be  applied  as  source  bias  for  these  regions.  The  fev 
station  biases  ve  have  obtained  for  sites  in  tectonic  regions  (Alaska,  Japan, 
South  America,  Nev  Guinea,  and  the  Caribbean)  are  hovever  mostly  positive 
(Figure  13)  and  none  are  less  than  -0.2  m^  units.  These  stations  are 
unfortunately,  as  mentioned  previously,  mostly  situated  relative  to  teleseismic 
activity  (i.e.,  A ^ 21^  as  defined  here)  such  that  the  knovn  zones  of  high 
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attenuation  in  their  vicinity  are  poorly  sampled.  Thus  we  cannot  estimate 
source  biases  with  any  degree  of  accuracy  (unless  we  accept  a relation  to 
velocities).  This  means  further  that  even  the  m^  values  as  corrected  for  bias 
do  not  necessarily  reflect  *true*  corrected  to  the  source,  which  can  only  be 
estimated  given  a knowledge  of  source  region  biases  (and  also  the  effects  of 
station  detection  characteristics). 

The  biases  we  have  obtained  from  large  events  (more  than  15  reports)  have 
been  applied  to  the  entire  set  of  events  (defined  such  that  at  least  3 of  the 
72  stations  chosen  reported  m^  values),  and  the  average  m^  values  recalculated. 
Figure  15  shows  the  magnitude  frequency  distribution  of  these  38,316  events 
prior  to  and  after  application  of  the  biases  as  station  corrections.  There  is 
little  change  in  the  distribution  of  events  for  m^  > 5*55  but  at  lower 
magnitudes  a substantial  redistribution  has  taken  place  towards  higher 
magnitudes.  This  is  caused  by  the  disproportionate  contribution  of  western 
U.S.  stations,  all  with  negative  biases,  to  the  total  number  of  station 
magnitude  reports,  and  will  be  discussed  in  the  next  section. 

We  now  wish  to  test  whether  the  application  of  these  biases  has  improved 
the  accuracy,  or  decreased  the  scatter,  of  the  individual  station  observations 
relative  to  the  average  for  a particular  event.  Clearly  it  will  do  so  for 
the  larger  events  (nearly  all  of  m^  > 5*0)  from  which  we  have  derived  the 
biases . 

We  define  the  ^scatter*  of  station  observations  relative  to  the  average 
by 
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mb 


Fig,  15.  Incremental  magnitude-frequency  curves  (as  defined  in  Fig.  3) 
for  38,316  events  over  196U-73  with  m^^  values  reported  by  3 or  more  of 
the  72  stations  chosen,  prior  to  (a)  and  after  (b)  application  of  mean 
biases  as  station  magnitude  corrections. 


for  events  in  a particular  magnitude  range  rri^  ■>  ra^  + where  is  the  total 

number  of  events  in  the  magnitude  range,  with  associated  station 

observations,  there  being  N.  (>3)  observations  m.  for  each  event  of  average 

^ ^ J 

magnitude  m^ . 

Figure  l6  shows  the  scatter  as  defined  above  for  Am^  = 0.1  rri^  unit  classes 
for  3.0  < average  rri^  < 6.7  for  both  the  uncorrected  and  corrected  datasets.  We 
can  see  that  application  of  the  biases  has  decreased  the  scatter  by  at  least 
15^  in  the  range  ^.3  - 6.0  m^.  Below  = ^-0  the  scatter  has  been  increased: 
this  is  presumably  because  the  effects  of  detection  characteristics  are 
particularly  severe  here.  Clearly  we  could  have  further  reduced  the  scatter  by 
applying  different  station  biases  for  each  seismic  region  (and  year!). 
Unfortunately  we  could  find  little  rational  explanation  for  the  yearly 
variation  in  station  bias,  and  no  consistency  in  regional  biases  across  small 
receiver  regions;  these  problems  should  be  solved  before  more  sophisticated 
bias  corrections  are  applied. 


U2 


SCATTER 


Fig.  l6.  Scatter,  as  defined  in  the  text,  as  a function  of  m-i^  (each  point 
denotes  scatter  for  all  events  of  mt)  + 0.1 ) before  (a),  and  after  (b) 

application  of  station  corrections.  Only  events  reported  by  3 or  more 
of  72  stations  are  used. 
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X.  STATION  DISTRIBUTION  AND  APPARENT  CHANGES  IN  SEISMICITY 


Table  I shows  that  for  the  first  6 years  of  the  time  period  196^-73 

considered  here  a very  large  proportion  of  all  station  reports  were  from 

stations  in  the  U.S.  This  proportion  declined  considerably  over  1970-73.  In 

particular,  5 Vela  stations  (BMO,  CPO,  TFO,  UBO,  and  WMO),  which  were  operated 

as  short  period  arrays,  and  two  other  stations  in  the  Western  U.S.  (DUG,  EUR) 

contributed  a vast  number  of  station  reports  in  earlier  years.  WMO  and  CPO, 

although  in  the  eastern  U.S.  province  of  Figure  6,  have  negative  biases  unlike 

other  stations  in  the  same  region:  this  may  be  partly  due  to  the  response 

characteristics  of  the  instruments,  which  are  capable  of  recording  much  shorter 

20 

period  signals  than  normal  (e.g.,  WWSSN ) stations  . Division  of  reported 

amplitude  A by  dominant  period  T does  not  entirely  compensate  for  the  much 

5 

higher  attenuation  at  shorter  periods.  Table  1 of  Evernden  and  Clark  also 
appears  to  show  that  magnitude  biases  for  WMO  and  CPO  are  anomalously  low 
compared  to  other  stations  in  the  Eastern  U.S.  province.  The  effect  of  these 
stations,  and  others  in  the  Western  U.S.  which  have  large  negative  biases,  is  a 
noticeable  reduction  in  average  in^  values,  x^^^'^icularly  at  lower  magnitudes. 
This  has  been  noted  in  the  previous  section  in  the  discussion  of  F’igure  15. 
Figure  17  shows  the  total  number  of  station  reports  per  year  for 

(i)  all  stations,  events  of  N^^  ^ 3,  N^^=station  reports/ 
event 

(ii)  72  stations,  events  of  N , >3 

st 

(iii)  ih  stations  in  continental  U.S.  (7  stations  below 
plus  ALQ,  BOZ,  GOL,  LAO,  LON,  NEW,  TUC) 


14)4 


STATION  nrih  REPORTS  (xlOOO) 


Fig.  17.  Total  station  m-^  reports/year  from 


(a)  all  stations 

(b)  72  stations  chosen 

(c)  US  stations 

(d)  5 Vela  stations  plus  DUG,  EUR. 


(iv)  the  7 stations  BMO,  CPO,  TFO,  UBO,  WMO,  DUG,  and 


EUR. 

During  196^4-69  stations  in  category  (iii)  above  contributed  from  to  6^^%  of 
all  station  observations,  and  those  in  category  (iv)  from  3^  to  58%.  These 
proportions  decline  to  ^>^30%  and  ^^^20%  respectively  in  1970-73.  Some  of  the  Vela 
stations  ceased  to  operate  as  arrays  (UBO,  CPO)  or  stopped  reporting  altogetlier 
(TFO,  WiO)  in  1969-70. 

The  reporting  performance  of  the  Vela  stations  is  truly  remarkable,  as 

shown  in  Figure  l8,  - BMO  never  reported  less  than  50%  of  all  events  of  N 3 

st 

for  our  72  station  data  set  in  any  one  given  year  and  in  I96U-65  UBO  reported 
over  80%  of  all  events.  The  number  of  events  reported  by  at  least  3 stations 
(all  stations,  not  just  the  72  chosen  here)  has  varied  (see  last  row  of  Table 
I)  from  3170  to  5132*,  the  lowest  number  is  for  1970,  when  some  Vela  stations 
stopped  or  reduced  their  reporting;  the  highest  is  for  1965,  in  which  the  Rat 
Island  sequence  contributed  over  1000  events. 

Figures  19  and  20  illustrate  the  variation  in  magnitude  distribution  due 
to  the  changing  contribution  with  time  of  the  7 stations  of  category  (iv)  above 
with  low  biases.  Figure  19a  shows  the  average  number  of  events /year, 
uncorrected  for  bias,  for  I96U-69  and  1970-73.  As  well  as  a considerable 
^‘eduction  in  the  number  of  smaller  events  for  the  later  period  (and  in  the 
annual  number  of  events  from  U096  to  3^33),  there  is  also  a redistribution  of 
events  towards  higher  magnitudes.  Application  of  the  biases  (Figure  19b) 
reduces  the  apparent  difference  in  seismicity  for  these  two  time  periods 
considerably,  though  there  is  still  a marked  reduction  in  the  number  of  smaller 
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PERCENTAGE  OF  TOTAL  GLOBAL  EVENTS  REPORTED 


|l8-2-  12938] 


Fig.  l8.  Percentage  of  total  global  events 
(3  or  more  of  72  stations  reporting  m^)  for 
which  values  have  been  given  by  5 Vela 
stations,  DUG,  and  EUR,  per  year. 
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Fig,  19.  Incremental  magnitude  frequency  curves,  reduced  to  annual  rates 
for  196^-69,  1970-73,  prior  to  (a'^  and  after  (b)  application  of  station 
m-K  corrections. 


mb  rOb 

Fig.  20.  Ratios  of  annual  seismicity,  I96U-69  and  1970-73,  to  annual 
seismicity,  196U-73,  as  a function  of  m-j^  (O.l  magnitude  classes)  be- 
fore and  after  station  corrections  have  been  applied. 


events  for  19T0-T3  compared  to  196^4-69.  Figure  20  expresses  the  difference  in 
another  manner,  and  shows  the  ratios  of  the  annual  number  of  events  in  a given 
0.1  range  for  6^4-69  and  70-73  to  that  for  the  entire  time  period  196^4-73. 

The  difference  from  mean  seismicity  (196^4-73)  for  in^  = ^4.8  is  reduced  by  a 
factor  of  2 on  application  of  the  biases  as  station  m^  corrections.  There  is 
still  an  apparent  increase  in  seismicity  for  the  later  time  period:  this  is 
possibly  because  the  negative  biases  for  Western  U.S.  stations  are 
underestimated  since  the  average  event  m^  for  I96U-69  is  reduced  by  the 
disproportionate  contribution  of  these  stations.  It  is  clear  that  network 
detection  capability  has  been  seriously  degraded  by  the  closure  or  reduction  in 
reporting  ability  of  these  Vela  stations.  It  is  somewhat  paradoxical  that 
although  there  are  more  stations  reporting  in  later  years  they  do  not  do  as 
well  at  lower  magnitudes  as  the  fewer  stations  of  196^4-69. 

Figure  21  gives  the  number  of  events  of  > 3.0,  ^4.0,  ^.5,  5*0,  5.5n  and 

6.0  reported  in  each  year.  It  can  be  seen  that  with  the  exception  of  a obvious 
21 

low  in  1966-67  the  number  of  events  of  m^  > 5*0  has  remained  effectively 
constant  as  measured  by  our  72-station  network,  and  m^  = 5 is  clearly  close  to 
the  detection  ability  of  the  network  for  all  years.  Note  that  the  aftershocks 
of  the  Alaskan  earthquake  of  196^4,  and  the  Rat  Island  sequence  of  1965,  have 
not  been  removed  from  the  seismicity.  Application  of  the  biases  has  not 
substantially  changed  the  annual  numbers  of  events  of  nr  > 5.0. 
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YEAR 


Fig.  21.  Annual  number  of  events,  reported  by  3 or  more  stations 
of  72-station  network,  m]^  > U.O,  4.5,  5.0,  5*5  and  6.0,  before 
and  after  station  corrections  have  been  used. 
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XI.  APPLICATIONS  TO  SEISMIC  DISCRIMINATION 


Use  of  the  station  biases  given  here  will  probably  reduce  the  scatter  in 
the  diagram,  though  the  biases  are  not  large  enough  to  greatly  improve 

the  separation  of  the  earthquake  and  explosion  populations.  If  one  also  takes 
into  account  biases  introduced  in  the  source  region  (which  we  have 
unfortunately  been  unable  to  determine  because  of  the  lack  of  stations  in 
seismic  regions)  the  improvement  could  be  dramatic. 

The  application  of  both  source  and  receiver  biases  can,  as  demonstrated  by 
Marshall^^,  drr:imatically  improve  the  linearity  of  the  magnitude-yield  relation. 
The  combined  effects  of  source  and  receiver  biases  is  such  that,  for  example, 
an  explosion  of  given  yield  detonated  in  the  Lake  Baikal  region  (rift  structure 
and  therefore  presumably  large  negative  bias)  and  recorded  in  the  Western  U.S., 
will  give  an  rri^  up  to  0.8  units  lower  than  that  for  one  of  the  same  yield  in 
Eastern  Kazakhstan  (stable  aseismic  platform)  reported  in  Canada  and  the 
Eastern  U.S.  Station  and  source  may  well,  in  certain  instances,  be 
geographically  located  such  that  such  large  variations  in  m^  for  a given  yield 
may  become  a reality. 

The  results  of  this  study  of  bias  also  demonstrate  the  dangers  of 
assigning  event  magnitudes  on  the  basis  of  one  station  report  only  (as  done, 
for  example,  by  the  NEIS  Earthquake  Data  Reports).  For  a station  whose  biases 
with  respect  to  the  mean  of  a large  network  of  stations  are  distributed 
normally  with  mean  -0.3  and  standard  deviation  0.35  m^  units  (typical  Western 
U.S.  values)  there  exists  a probability  of  0.50  that  the  station  magnitude  will 
be  lower  by  0.3  m^  units,  and  a probability  0.20  that  it  will  be  lower  by  0.6, 


than  the  average  of  a large  network  observing  the  same  event.  This  is  one  of 

the  justifications  for  studying  here  only  events  for  which  has  been  reported 

by  3 or  more  stations  (except  that  if  all  these  3 are  in  the  Western  U.3.  we 

are  no  wiser!).  A more  comprehensive  statistical  analysis  of  taking  into 

account  random  variations,  biases,  and  the  effects  of  detection  thresholds, 

19  11 

such  as  that  proposed  by  Ringdal  and  Christoff ersson  et  al.  is  clearly 
desirable . 
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XII.  CONCLUSIONS 


We  have  shown,  using  a network  of  the  best  12  stations  for  196^-73,  that 
there  exist  substantial  station  biases  (up  to  at  least  0.^  units)  as 
measured  relative  to  average  iri^.  These  biases  are  well  correlated  with 
tectonic  structure  and  previous  measurements  of  lateral  variations  in 
attenuation.  There  is  also  some  evidence  that  bias  is  correlated  with  P- 
velocity  in  the  upper  mantle  Application  of  these  biases  to  the 

calculation  of  average  la^  for  a particular  event  reduces  the  scatter  in  the 
average  and  also  removes  many  of  the  apparent  changes  in  seismicity  with 
time,  despite  the  fact  that  the  biases  themselves  are  clearly  a function  of 
source  region  (and  occasionally  time!).  If  the  seismic  research  observatory 
(SRO)  stations  now  being  installed,  which  have  much  the  sajne  fairly  broad-band 
short  period  response  as  the  Vela  stations,  are  operated  and  read  as  well  as 
the  latter,  their  detection  capability  should  be  much  better  than  that  of  the 


existing  global  station  network. 
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APPENDIX  A 


GEOGRAPHICAL  LOCATIONS  OF  STATIONS  USED  IN  TABLES  II  AND  IV 


Asterisk  (*)  denotes  array  operation  for  some  or  all  of  196^1-73 ; 
(W)  denotes  WWSSN  station. 


ABU 

Abuyama,  Honshu,  Japan 

ALE 

Alert,  Northwest  Territories,  Canada 

ALQ 

(w) 

Albuquerque,  New  Mexico 

ASP 

Alice  Springs,  Northern  Territory,  Australia 

ATL 

(w) 

Atlanta,  Georgia 

BHA 

Broken  Hill,  Zambia 

BKS 

(w) 

Berkeley,  California 

BMO* 

Blue  Mountains,  Oregon 

BNG 

Bangui,  Central  African  Republic 

BNS 

Ben sb erg.  West  Germany 

BOZ 

Bozeman,  Montana 

BRG 

Berggiesshubel,  East  Germany 

BUH 

Buhlerhohe , West  Germany 

BUL 

(w) 

Bulawayo,  Rhodesia 

CAN 

Canberra,  Capital  Territory,  Australia 

CAR 

(w) 

Caracas,  Venezuela 

CIR 

Chiredzi , Rhodesia 

CLL 

Collmberg , East  Germany 
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CLK 

Chileka,  Malawi 

COL  (W) 

College,  Alaska 

COP  (W) 

Copenhagen , Denmark 

CPO* 

Cumberland  Plateau,  Tennessee 

CPO* 

Cumberland  Plateau,  Tennessee 

DAG  (W) 

Danmarks  Havn,  Greenland 

DUG  (W) 

Dugway,  Utah 

EDM 

Edmonton,  Alberta,  Canada 

EKA* 

Eskdalemuir , Scotland 

ESK  (W) 

Eskdalemuir , Scotland 

EUR 

Eureka,  Nevada 

EEC 

FI in  Flon,  Manitoba,  Canada 

ELO  (W) 

Florissant,  Missouri 

EUR 

Furstenfeldbruck,  West  Germany 

GBA* 

Gauribidanur , India 

GDH  (W) 

Godhavn , Greenland 

GIL 

Gilmore  Creek,  Alaska 

GOL  (W) 

Golden,  Colorado 

GRE* 

Grafenberg,  West  Germany 

HES* 

Hagfors , Sweden 

HYB 

Hyderabad , India 

I LG 

Inge  Lehmann,  Greenland 

INK 

Inuvik,  Northwest  Territories,  Canada 

KBL  (W) 

Kabul , Afghanistan 

57 


KEV  (W) 

Kevo , Finland 

KHC 

Kasperske  Hory,  Czechoslovakia 

KJF 

Kajaani,  Finland 

KJN 

Ka j aani , Finland 

KLG 

Kalgoorlie,  Western  Australia 

KOD  (W) 

Kodaikanal,  India 

KON  (W) 

Kongsberg , Norway 

KRA 

Krakow,  Poland 

KRR 

Karoi,  Rhodesia 

KTG  (W) 

Kap  Tobin,  Greenland 

LAO* 

LASA,  Montana 

LHN 

Lillehammer , Norway 

LJU 

L j ubl j ana , Yugoslavia 

LON  (W) 

Longmire , Washington 

LOR  (W) 

Lormes,  France 

LPB  (W) 

La  Paz,  Bolivia 

LPS  (W) 

La  Palma,  El  Salvador 

MAW 

Maw son , Antarctica 

MBC 

Mould  Bay,  Northwest  Territories,  Canada 

MOX 

Moxa,  East  Germany 

MTD 

Mount  Darwin,  Rhodesia 

MUN  (W) 

Mundaring,  Western  Australia 

NAO* 

NORSAR , Norway 

NDI  (W) 

New  Delhi,  India 
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NEW  (W) 

Newport,  Washington 

NIE 

Niedzica,  Poland 

NOR  (W) 

Nord,  Greenland 

NP- 

North  Pole,  Northwest  Territories,  Canada 

NUR  (W) 

Nurmi jarvi , Finland 

OIS 

Oishiyama,  Honshu,  Japan 

PAE 

Paea,  French  Polynesia 

PMG  (W) 

Port  Moresby,  Papua 

PMO 

Pomariorio,  French  Polynesia 

PMR 

Palmer , Alaska 

PNS  (W) 

Penas,  Bolivia 

PNT 

Penticton,  British  Columbia,  Canada 

POO  (W) 

Poona,  India 

PPN 

Papenoo , French  Polynesia 

PPT 

Papeete,  French  Polynesia 

PRE  (W) 

Pretoria,  South  Africa 

PRU 

Pruhonice,  Czechoslovakia 

RAB  (W) 

Rabaul , New  Britain 

ROD  (W) 

Rapid  City,  South  Dakota 

RES 

Resolute  Bay,  Northwest  Territories,  Canada 

RIV  (W) 

Riverview,  New  South  Wales,  Australia 

ROL 

Rolla,  Missouri 

RUV 

Rauvai,  French  Polynesia 

TFO* 

Tonto  Forest,  Arizona 
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TPT 

Tiputa,  French  Polynesia 

TRN 

(w) 

Trinidad,  Trinidad  and  Tobago 

TSK 

Tsukuba,  Honshu,  Japan 

TUC 

(W) 

Tucson,  Arizona 

TUL 

(W) 

Tulsa,  Oklahoma 

TVO 

Taravao,  French  Polynesia 

UBO* 

Uinta  Basin,  Utah 

VAH 

Vaihoa,  French  Polynesia 

WIN 

(w) 

Windhoek,  Namibia 

WMO* 

Wichita  Mts.  Oklahoma 
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